Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex 
Avoidance System (Wake VAS) Prediction Algorithm 

NAG 1-03043 Final Research Grant Report 
NASA-Langley Research Center 
Hampton, Virginia 23681-2199 
ATTN: Fred H. Proctor 


Michael L. Kaplan* and Yuh-Lang Lin 
Department of Marine, Earth and Atmospheric Sciences 
North Carolina State University 
Raleigh, North Carolina 27695-8208 


December 2004 


*Michael L. Kaplan, Department of Marine, Earth and Atmospheric Sciences, North 
Carolina State University, Raleigh, North Carolina 27695-8208, E-MAIL: 


1 



mlkaplan@unity.ncsu.edu , PHONE:919-515-1442, FAX:919-48 1-9332. 

1. Summary of the Research 

During the research project, sounding datasets were generated for the region 
surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, 
Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los 
Angeles, CA. The numerical simulation of winter and summer environments during which 
no instrument flight rule impact was occurring at these 9 terminals was performed using 
the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) 
model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical 
resolution in the planetary boundary layer (e.g, note Appendix A). The soundings from the 
1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical 
dependent variables as well as derived quantities, i.e., 3-dimensional wind components, 
temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation 
rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above 
ground level. After partial validation against field experiment datasets for Dallas as well as 
larger scale and much coarser resolution observations at the other 8 airports, these 
sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling 
program. The application of these datasets being to determine representative airport 
weather environments to diagnose the response of simulated wake vortices to realistic 
atmospheric environments. These virtual datasets are based on large scale observed 
atmospheric initial conditions that are dynamically interpolated in space and time 
employing a numerical simulation model based on a comprehensive hydro- 
thermodynamical equation set. The 1 km nested-grid simulated datasets providing a very 
coarse and highly smoothed representation of airport environment meteorological 
conditions. Details concerning the airport surface forcing are virtually absent from these 
simulated datasets although the observed background atmospheric processes have been 
compared to the simulated fields and the fields were found to accurately replicate the 
flows surrounding the airport where coarse verification data were available as well as 
where airport scale datasets where available at Dallas (e,g., note Appendix A). The 
simulated profiles of atmospheric dependent variables can be viewed online at the 
following links: 

1 ) http://mesolab.meas.ncsu.edu/%7Eamhoggarth/wakevortex/wake research.html , 2) 
http://mesolab.meas.ncsu.edu/~mtk3581/wake/ , 3) http://mesolab.meas.ncsu.edu/~paul/ 


2 


Appendix A 


Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex 
Advisory System - Dallas-Fort Worth Wake-VAS Preliminary Simulation Experiments 


NAG1-0343 Research Task Report for Summer 2004 
NASA-Langley Research Center 
Hampton, Virginia 23681-2199 
ATTN: David K. Rutishauser 


Michael L. Kaplan* and Yuh-Lang Lin 
Department of Marine, Earth and Atmospheric Sciences 
North Carolina State University 
Raleigh, North Carolina 27695-8208 


September 2004 


3 



*Michael Kaplan, Department of Marine, Earth and Atmospheric Sciences, North Carolina 
State University, Raleigh, NC, 27695-8208, E-MAIL:mlkaplan@unity.ncsu.edu, 
PHONE:919-5 15-1442, FAX:919-515-1683. 

1. Task Overview 

This report summarizes research on the development of virtual atmosphere sounding 
data sets for the Wake-Vortex Advisory System (Wake-VAS) during the summer of 2004. 
During this period the focus was on generating virtual sounding data sets at Dallas-Fort 
Worth International Airport (DFW) during the September 1997 AVOSS deployment. 
Soundings generated by the NHMASS 6.2 model, i.e., the Terminal Area PBL Prediction 
System (TAPPS), where directly compared to MIT Lincoln Laboratory synthesized 
observed soundings of u and v wind components as well as observations of 40 m eddy 
dissipation rate and turbulence kinetic energy at DFW. The comparison would represent 
evaluations of potential control simulations in anticipation of more complex future 
simulation studies. This represents precursor analyses in preparation for future simulation 
experiments in which local observations would be assimilated into the predictive system. 
These simulated soundings have highly detailed vertical structure for use in diagnosing 
realistic vertical wind shears, temperature lapse rates and, therefore, potential turbulence 
generating features in the planetary boundary layer at the DFW airport. The simulations 
were validated against the MIT observations to diagnose the status of their existing utility 
as components of a future Wake-VAS predictive system. 

2. Modeling System/Numerical Simulation Experiments 

The numerical model, NHMASS 6.2, is a nonhydrostatic and terrain following 
simulation model. This is a nonhydrostatic version of the Terminal Area PBL Prediction 
System (TAPPS) (Kaplan et al. 2000). The model lid was set at 100 hPa with focused 
vertical resolution below 1000 m elevation above ground level (AGL). The version of 
NHMASS employed for these simulation experiments included a 1.5 order prognostic 
turbulence kinetic energy planetary boundary layer formulation after Therry and Lacarrere 
(1983) and a convective parameterization scheme after Kain and Fritsch (1992). A 
comprehensive surface energy budget and soil moisture predictive scheme is utilized in the 
numerical model so that diurnal radiative forcing of the earth’s surface and hydrological 
forcing is included in the surface layer of the simulated atmosphere (Mahrt and Pan 1984). 
Surface soil moisture was set at a climatological value for September in Dallas and land 
use/land cover databases as well as surface albedo where also generated based on 
prescribed datasets (Anderson et al. 1976; Noilan and Planton 1989). The eddy dissipation 
rate formulation (e) employed in diagnostic calculations is consistent with that utilized in 
the turbulence kinetic energy boundary layer formulation. It is represented by the 
following expression from Bougeault and Lacarrere (1989) and Therry and Lacarrere 
(1983): 


£ = Ce 3/2 /Xe 


( 1 ) 


where: 
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e = first guess turbulence kinetic energy derived from the full time dependent equation 


O = a first order numerical coefficient 


• • = a characteristic length scale of the energy-containing eddies 


The value of this characteristic length scale (• •) is primarily a function of Richardson 
number wherein numerous thresholds exist for differing stability and shear regimes in 
order to differentiate the mixing length of eddies in stable, convective and transitional 
boundary layers. 

The numerical experiments consisted of performing several one-way doubly-nested 
numerical simulations with NHMASS 6.2 during the September 1997 AYOSS deployment 
in the Dallas-Fort Worth, Texas metropolitan region. The simulations were centered on 
DFW. Three numerical grids were employed for these experiments: 1) 36 km horizontal 
resolution, 2) 6 km horizontal resolution and 3) 1 km horizontal resolution. Matrices 
included a 100 x 100 x 59, 50 x 50 x 59 and 25 x 25 x 59 grid, respectively, with very 
detailed and identical vertical resolution within the first 1000 m AGL of the numerical 
model for all 3 grids. The time periods of integration for the 3 grids were 24 hours, 24 
hours and 30 minutes regenerated for the entire 24-hour period, respectively. This means 
that every 30 minutes a 1 km horizontal resolution simulation was performed from the 6 
km simulated initial and lateral boundary condition data. The 6 km initial and lateral 
boundary conditions were derived from the 36 km simulation, which, in turn received its 
initial and lateral boundary conditions from NWS ETA gridded analyses, which were 
reanalyzed with synoptic scale rawinsondes as well as ASOS surface observations. No 
local (special AVOSS) DFW observations were utilized in the numerical simulations. The 
initial and lateral boundary condition datasets were solely derived from the NWS ETA 
analyses gridded fields as well as reanalyzed standard rawinsonde and aviation surface 
datasets. The simulations were initialized at 8 different time periods, i.e., 0000 UTC 16, 

17, 18 and 19 September 1997 (group #1) and 1200 UTC 16, 17, 18 and 19 September 
1997 (group #2). The groupings would allow statistical comparisons of the relative 
accuracy of the simulations comparing 0000 and 1200 UTC initialization times. Vertical 
soundings were constructed centered on DFW every 30 minutes from the 1 km 
simulations and interpolated to 5 m vertical resolution for comparison to the available 
observed fields at 200 levels up to 1000 m AGL. These simulated soundings included the 
following variables: 1) pressure (hPa), 2) height (m AGL), 3) temperature (K), 4) dew 
point temperature (K), 5) u wind component (m/s), 6) v wind component (m/s), 7) w wind 
component (m/s), 8) turbulence kinetic energy (m 2 /s 2 ) and 9) eddy dissipation rate (m 2 /s 3 ). 
These were sent to NAS A-Langley for potential use in the Virtual Air Space Modeling 
Program. Figure 1, below, represents an example of one of these 384 soundings. 


P (mb) Z (m) T (K) TD (K) U(m/s) V(m/s) W(m/s) TKE TKED 


995.142 

5.00000 

300.615 

294.579 

-2.29894 

-.389851 

— . 1369' 7E-01 

0 . 100000E-03 

0.00000 

994.577 

10.0000 

300.622 

294.568 

-2.30653 

-.332376 

- . 144184E-01 

0.281843E-01 

0 . 431458E-03 

994.012 

15.0000 

300.631 

294.551 

-2.30987 

-.240490 

- . 152 993E-01 

0 . 574212E-01 

0 . 849109E-03 

993.447 

20.0000 

300.644 

294.526 

-2.30443 

-.105129 

- . 16157 9E-01 

0.758775E-01 

0 . 102403E-02 

992.883 

25.0000 

300.661 

294.491 

-2.28772 

0 . 7 61115E-01 

- . 169777E-01 

0.8217 31E-Q1 

0 . 934865E-03 

992.319 

30.0000 

300.672 

294.454 

-2.25700 

0.269599 

- . 177050E-01 

0 . 809197E-01 

0 . 731693E-03 

991.755 

35.0000 

300.666 

294.421 

-2.20826 

0.442406 

- . 182775E-01 

0 . 7 6187 0E-01 

0 . 555970E-03 
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991.192 

40.0000 

300.646 

294.394 

-2.14447 

0.587835 

- . 187137E-01 

0.711182F-01 

0 . 448001E-03 

990.629 

45.0000 

300.617 

294.370 

-2.07236 

0.713316 

- . 1907 68E-01 

0 . 685720E-01 

0 . 392788E-03 

990.066 

50.0000 

300.585 

294.350 

-1.9971.7 

0.823573 

- . 194162E-01 

0.714182E-01 

0 . 381731E-03 

989.504 

55.0000 

300.550 

294.331 

-1 . 92102 

0.923215 

- . 197466E-01 

0 . 7 94071E-01 

0 . 402202E-03 

988.942 

60.0000 

300.512 

294.314 

-1.84439 

1.01648 

- . 200634E-01 

0.908216E-01 

0 . 440239E-03 

988.380 

65.0000 

300.472 

294.297 

-1.76767 

1.10658 

- . 203600E-01 

0.104461 

0 . 484739E-03 

987.818 

70.0000 

300.428 

294.282 

-1.69116 

1.195'. 5 

- . 206360E-01 

0.119178 

0 . 529359E-03 

987.256 

75.0000 

300.383 

294.267 

-1.61515 

1.28225 

- . 208995E-01 

0.133612 

0 . 572433E-03 

986.695 

80.0000 

300.337 

294.252 

-1.53997 

1.36765 

- . 211575E-01 

0.146335 

0 . 612339E-03 

986.134 

85.0000 

300.293 

294.237 

-1.46603 

1.45094 

-.214152E-01 

0.156033 

0 . 647493E-03 

985.573 

90.0000 

300.249 

294.222 

-1.39333 

1.53228 

- . 216716E-01 

0.163891 

0 . 67 9738E-03 

985.013 

95.0000 

300.205 

294.207 

-1.32178 

1.61194 

- . 21 922 9E-01 

0.172154 

0.712358E-03 

984.452 

100.000 

300.158 

294.192 

-1.25147 

1.68989 

- . 221645E-01 

0.183138 

0 . 748705E-03 

983.892 

105.000 

300.107 

294.178 

-1 .18261 

1.76597 

- . 223914E-01 

0.198884 

0 . 7 91624E-03 

983.332 

110.000 

300.053 

294.164 

-1.11498 

1.84049 

- . 22 6055E-01 

0.217562 

0 . 837848E-03 

982.772 

115.000 

299.999 

294.149 

-1.04831 

1.91379 

- . 228103E-01 

0.235995 

0 . 881827E-03 

982.213 

120.000 

299.946 

294.135 

-.982475 

1.98602 

- . 230077E-01 

0.250960 

0 . 917 638E-03 

981.654 

125.000 

299.895 

294.121 

-.917500 

2.05719 

- . 231985E-01 

0.258848 

0 . 938368E-03 

981.095 

130.000 

299.847 

294.107 

-.853308 

2.12737 

- . 233827E-01 

0.259168 

0 . 942715E-03 

980.536 

135.000 

299.801 

294.092 

-.789728 

2.19679 

- . 235597E-01 

0.254631 

0 . 935077E-03 

979.977 

140.000 

299.755 

294.078 

-.726673 

2.26554 

- . 23727 9E-01 

0.247366 

0 . 918305E-03 

979.419 

145.000 

299.711 

294.063 

-.664144 

2.33360 

- . 238840E-01 

0.239060 

0 . 893774E-03 

978.861 

150.000 

299.665 

294.049 

-.602210 

2.40090 

- . 240237E-01 

0.231152 

0 . 8617 95E-03 

978.303 

155.000 

299.619 

294.035 

-.540820 

2.46747 

- . 241471E-01 

0.224016 

0 . 823924E-03 

977.745 

160.000 

299.573 

294.020 

-.479858 

2.53348 

- . 242571E-01 

0.217366 

0 . 782485E-03 

977.188 

165.000 

299.526 

294.006 

-.419275 

2.59897 

- . 243545E-01 

0.210973 

0 . 739267E-03 

976.631 

170.000 

299.480 

293.991 

-.359085 

2.66390 

- . 244382E-01 

0.204652 

0 . 6957 68E-03 

976.074 

175.000 

299.433 

293.977 

-.299358 

2.72818 

- . 245048E-01 

0.198232 

0.653421E-03 

975.517 

180.000 

299.386 

293.962 

-.240105 

2.79177 

- . 245537E-01 

0.191590 

0 . 612834E-03 

974 . 961 

185.000 

299.339 

293.946 

-.181210 

2.85483 

- . 245892E-01 

0.184664 

0 . 573457E-03 

974.405 

190.000 

299.293 

293.931 

-.122605 

2 . 91745 

- . 246140E-01 

0.177352 

0 . 534827E-03 

973.849 

195.000 

299.246 

293.916 

-. 642646E-01 

2.97965 

- . 246291E-01 

0.169488 

0 . 4 96603E-03 

973.293 

200.000 

299.200 

293.901 

-. 620795E-02 

3.04137 

- . 246338E-Q1 

0.160822 

0 . 458509E-03 

972.737 

205.000 

299.154 

293.886 

0.515062E-01 

3.10251 

- . 246262E-01 

0.151004 

0 . 420317E-03 

972.182 

210.000 

299.109 

293.872 

0.108893 

3.16308 

- . 246064E-01 

0.140254 

0 . 382435E-03 

971 . 627 

215.000 

299.063 

293.858 

0.166037 

3.22321 

- . 2457 68E-01 

0.129356 

0 . 345937E-03 

971.072 

220.000 

299.018 

293.843 

0.222990 

3.28298 

- . 245383E-01 

0 .119101 

0.312143E-03 

970.517 

225.000 

298.972 

293.829 

0.279772 

3.34241 

- . 244905E-01 

0.110384 

0 , 282812E-03 

969.963 

230.000 

298.925 

293.815 

0.336376 

3.40146 

- . 244313E-01 

0.104346 

0 . 2 60423E-03 

969.409 

235.000 

298.878 

293.801 

0.392764 

3.46003 

- . 243569E-01 

0.102560 

0 . 248570E-03 

968.855 

240.000 

298.830 

293.787 

0.448918 

3.51806 

- . 242643E-01 

0.106128 

0 . 249807E-03 

968.301 

245.000 

298.780 

293.772 

0.504933 

3.57568 

- . 241585E-01 

0,113385 

0 . 260190E-03 

967.748 

250.000 

298.731 

293.758 

0.560917 

3.63303 

- . 240442E-01 

0.122560 

0 . 275591E-03 

967.195 

255.000 

298.681 

293.744 

0.616970 

3.69024 

- . 239248E-01 

0.132190 

0 . 2 92 692E-03 

966.641 

260.000 

298.632 

293.729 

0.673204 

3.74736 

- . 238031E-01 

0.140889 

0 . 308451E-03 

966.089 

265.000 

298.584 

293.714 

0.729760 

3.80446 

- . 236811F.-01 

0.147116 

0 . 319607E-03 

965.536 

270.000 

298.536 

293.699 

0.786831 

3.86160 

- . 235610E-01 

0.148937 

0.322159E-03 

964.984 

275.000 

298.490 

293.683 

0.844569 

3.91883 

- . 234445E-01 

0.144869 

0 . 313053E-03 

964.432 

280.000 

298.446 

293.667 

0.902816 

3.97611 

- . 2332 99E-Q1 

0.136401 

0 . 295207E-03 

963.880 

285.000 

298.402 

293.651 

0.961409 

4.03340 

- . 232157E-01 

0,125123 

0 . 271752E-03 

963.328 

290.000 

298.359 

293.634 

1.02024 

4.09065 

- . 231000E-01 

0.112170 

0 . 244917E-03 

962.777 

295.000 

298.316 

293.617 

1.07925 

4.14781 

- . 229810E-01 

0 . 984097E-01 

0 . 21 6383F.-03 

962.226 

300.000 

298.273 

293.601 

1.13838 

4.20482 

- . 228562E-01 

0.845768E-01 

0 . 187558E-03 

961.675 

305.000 

298.230 

293.585 

1.19763 

4.26160 

- . 227220E-01 

0.713958E-01 

0 . 159804E-03 

961.124 

310.000 

298.187 

293.569 

1.25699 

4.31801 

- . 225737E-01 

0.596892E-01 

0 . 134 64 F.-03 

960.573 

315.000 

298.144 

293.553 

1.31644 

4.37398 

- . 224092E-01 

0 . 4 987 93E-01 

0 . 112862F.-03 

960.023 

320.000 

298.100 

293.538 

1.37587 

4.42953 

- . 222328E-01 

0.415432E-01 

0 . 937316E-04 

959.473 

325.000 

298.056 

293.522 

1.43518 

4.48466 

- . 22047 9E-01 

0 . 3437 68E-01 

0 . 7 6727 9E-04 

958.923 

330.000 

298.012 

293.508 

1.49424 

4.53927 

218567E-01 

0.282644E-01 

0 . 616725E-04 

958.374 

335.000 

297.968 

293.493 

1.55289 

4.59320 

- . 21 6601F.-01 

0 . 232 604E-01 

0 . 487007E-04 

957.825 

340.000 

297.923 

293.479 

1.61089 

4 . 64618 

- . 214 583E-01 

0 . 195981E-01 

0 . 382840E-04 

957.276 

345.000 

297.878 

293.465 

1.66790 

4.69775 

- . 212506F.-01 

0.177269E-01 

0.313063E-04 

956.727 

350.000 

297.831 

293.451 

1.72343 

4.74728 

-, 21035 5E-01 

0 . 183854E-01 

0 . 292058E-04 

956.178 

355.000 

297.783 

293.439 

1.77683 

4.79389 

- . 208105E-01 

0.225683E-01 

0 . 339155E-04 

955.630 

360.000 

297.733 

293.427 

1.82832 

4.83784 

- . 205772E-01 

0 . 2 99181E-01 

0.447712F-04 

955.082 

365.000 

297.683 

293.416 

1.87860 

4.88004 

- . 203397E-01 

0 . 392834E-01 

0 . 595899E-04 

954.534 

370.000 

297 . 632 

293.405 

1.92820 

4 . 92112. 

-.201013E-01 

0 . 4 98135E-01 

0 . 7 67 803E-04 

953.986 

375.000 

297.581 

293.394 

1.97754 

4.96156 

- . 198652E-01 

0 . 608451E-01 

0 . 951335E-04 

953.439 

380.000 

297.531 

293.382 

2.02701 

5.00169 

- . 196353E-01 

0 . 717 963E-01 

0 . 113628F.-03 

952.891 

385.000 

297.481 

293.369 

2.07700 

5.04184 

- . 194160E-01 

0.820719E-01 

0 . 1312 67E-03 

952.344 

390.000 

297.433 

293.356 

2.12799 

5.08231 

- . 192137F.-01 

0.909679E-01 

0 . 14 6911E-03 

951.797 

395.000 

297.386 

293.340 

2.18062 

5.12348 

- . 19037 6E-01 

0.975575E-01 

0 . 159105E-03 

951.251 

400.000 

297.342 

293.321 

2.23582 

5.16580 

- . 189009E-01 

0.100545 

0.165836E-03 

950.704 

405.000 

297.302 

293.298 

2.29414 

5.20959 

- . 188116E-01 

0.991236E-01 

0 . 165849E-03 

950.158 

410.000 

297.265 

293.273 

2.35475 

5.25450 

- . 187561E-01 

0 . 944914E-01 

0 . 160952E-03 

949.613 

415.000 

297.230 

293.245 

2.41690 

5.30022 

- . 187216E-01 

0 . 877460E-01 

0 . 152801E-03 

949.067 

420.000 

297.196 

293.216 

2.48008 

5.34656 

- . 186987E-01 

0.796446E-01 

0 . 142525E-03 

948.522 

425.000 

297.164 

293.185 

2.54394 

5.39343 

- . 1867 98E-01 

0.707253E-01 

0 . 130912E-03 

947 . 976 

430.000 

297.134 

293.151 

2.60826 

5.44085 

- . 186582E-01 

0 . 613944E-01 

0..118535E-03 

947.432 

435.000 

297.106 

293.116 

2 . 67287 

5.48892 

- . 186264E-01 

0.519920E-01 

0 . 105841E-03 

946.887 

440.000 

297.082 

293.076 

2.73767 

5.53781 

- . 185755E-01 

0 . 428465E-01 

0 . 932285E-04 

946.343 

445.000 

297.062 

293.032 

2.80257 

5.58787 

- . 184 935E-01 

0 . 34327 6E-01 

0 . 811081E-04 

945.799 

450.000 

297.049 

292.980 

2.86750 

5.63962 

- . 183634E-01 

0.269063E-01 

0 . 699698E-04 

945.255 

455.000 

297.046 

292 . 918 

2 . 93242 

5.69372 

- . 181642E-01 

0. 21134 6E-01 

0 . 603702E-04 

944.712 

460.000 

297.050 

292.846 

2 . 99726 

5.74987 

- . 179032F.-01 

0.168534E-01 

0.521773E-04 

944.169 

465.000 

297.061 

292.769 

3.06192 

5.80735 

-.176006E-01 

0 . 1357 99E-01 

0 . 449457E-04 
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Figure 1 - Simulated 5 -meter vertical resolution DFW virtual soundings of pressure 
(hPa), height AGL (m), temperature (K), dew point temperature (K), u wind component 
(m/s), v wind component (m/s), w wind component (m/s), turbulence kinetic energy (m 2 
/s 2 ), and eddy dissipation rate (m 2 /s 3 ) valid at 1200 UTC 18 September 1997. 


3. Simulation Validation 


The validation datasets include DFW synthesized wind soundings, which represented a 
synthesis of DFW wind profiler, lidar, sodar and anemometer observations, as well as 5 m 
and 40 m observed eddy dissipation rate estimates prepared by the MIT Lincoln 
Laboratory (e.g., Dasey et al. 1998). The validation procedure included calculating both 
absolute errors and root mean square errors at times and levels where direct validation 
data from MIT where available. This required interpolating MIT Lincoln Laboratory 
sounding datasets to the same vertical resolution as the TAPPS simulated virtual 
soundings. 

Figures 2-13, below, depict tables of averaged absolute errors (absolute values of errors) 
of u and v wind component up to 600 m AGL, eddy dissipation rate at 40m AGL and 
turbulence kinetic energy at 40 m AGL over each group of 4 case studies, i.e., the 1200 
UTC and 0000 UTC initialized simulations. Only two time periods of observed data were 
employed to initialize these simulations, which are validated out through 24 hours. These 
errors are for a very small sample of case studies and represent a point verification for a 
group of similar weakly-forced large scale environments. Additionally, the MIT profiles 
may have systematic biases as they represent a synthesis of numerous sensor-derived 
observations with their own intrinsic biases and errors. The u and v wind component 
versus height validation (Figures 2-5) represent direct comparison of averaged 
observations versus averaged simulated fields over each group of 4 case studies at specific 
heights for all times. The u and v wind component versus time validation (Figures 6-9) 
represent averaged observations versus averaged simulated fields over each group of 4 
case studies for all heights at specific times. Note that the observed datasets were missing 
at times specified on each diagram. Only 40 m turbulence kinetic energy and eddy 
dissipation rate errors could be calculated since the 5 m simulated fields were not available 
due to their being below the first model level of ~7 m AGL. 

Specific findings of these small numbers of case studies can be summarized by the 
following general conclusions concerning u and v wind component errors versus height 
and time: 

1) The u component errors exhibited little height dependence while the v component 
errors exhibited significant relative height dependence from both groups. 
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2) The 1200 UTC initialized u and v component errors were much smaller than the 
0000 UTC initialized u and v component errors versus height and time. 

3) The u and v component errors are similar in magnitude when initialized at 1200 
UTC but the v component is much larger than the u component wind error 
initialized at 0000 UTC. 

4) The maximum error with respect to height for both groups and both components is 
at -400 m +- -100 m AGL. 

5) The maximum error with respect to time for both groups is during the 
midafternoon period. 

6) Both the u and v wind components exhibited double temporal error maxima from 
the 0000 UTC initialized simulations but only a single temporal maxima from the 
1200 UTC simulations. 

7) The u and v components exhibited a double temporal minima in the 1200 UTC 
simulations but not in the 0000 UTC simulations where a single temporal minima 
occurred. 

As far as eddy dissipation rate and turbulence kinetic energy errors are concerned (note 
Figures 10-13): 

1) Eddy dissipation rate errors peak ~10‘ 2 m 2 /s 3 after the first 30 minute period but 
then drop to ~10‘ 3 m 2 /s 3 for the remainder of the simulation for both groups. 

2) Turbulence kinetic energy errors are consistent with wind component errors in 
terms of their temporal variation. 

The most important general and preliminary findings to be derived from these 
experiments are the following: 

1) 1200 UTC initial data enabled a much more accurate simulation of noctumally- 
forced boundary layer jets as well as, possibly, the background synoptic flow 
regime at DFW. This is likely because the 1200 UTC large scale analyses 
contained more of a signal of these boundary layer features, in particular, which 
had not developed at 0000 UTC and where, therefore absent from the larger scale 
initial datasets at 0000 UTC. May also be, in part, a result of the times of the 
availability of observations. 

2) The v wind component variance is much more pronounced in the lower 
troposphere for nocturnal boundary layer jets than is the u wind component 
variance. 

3) The simulation of stable boundary layer may lend itself to deterministic simulation 
better than the convective boundary layer does as errors tended to peak during the 
period of most significant surface sensible heating and convective overturning 
within the boundary layer. 

4) Eddy dissipation rate simulation is likely very sensitive to the initial data and initial 
data adjustment processes employed in a deterministic numerical model as these 
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U Component error (m/s) 


errors peaked right after model initialization. Overall these errors were very large 
but are likely highly sensitive to the formulation employed for both the 
observations and simulations reflecting differences in these formulations. 


5) The most frequent level of the nocturnal boundary level jets at DFW is -400 m 
AGL. 



Figure 2 - Simulated averaged u wind component errors (m/s) versus height (m AGL) 
from the 1200 UTC initial data. 
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U-Component error vs. Height 



Height (m) 

Figure 3 - Simulated averaged u wind component error (m/s) versus height (m AGL) from 
the 0000 UTC initial data. 


V-Component error vs. Height 
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V-Component error (m/s) 


Figure 4 - Simulated averaged v wind component error (m/s) versus height (m AGL) from 
the 1200 UTC initial data. 


V-Component error vs. Height 



Figure 5 - Simulated averaged v wind component error (m/s) versus height (m AGL) from 
the 0000 UTC initial data. 
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U-Component error vs. Time 



Local DFW Time 


Figure 6 - Simulated averaged u wind component error (m/s) versus time (hr) from the 
1200 UTC initial data. 


U-Component error vs. time 



Local DFW Time 
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Figure 7 - Simulated averaged u wind component error (m/s) versus time (hr) from the 
0000 UTC initial data. 


V-Component error vs. Time 



Local DFW Time 


Figure 8 - Simulated averaged v wind component error (m/s) versus time (hr) from the 
1200 UTC initial data. 
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V-Component error vs. time 
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Figure 9 - Simulated averaged v wind component error (m/s) versus time (hr) from the 
0000 UTC initial data. 


EDR Error vs. Time 
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EDR error 


Figure 10 - Simulated averaged eddy dissipation rate error (m 2 /s 3 ) versus time (hr) at 40 
m AGL from the 1200 UTC initial data. 


EDR Error vs. Time 



Local DFW Time 


Figure 1 1 - Simulated averaged eddy dissipation rate error (m 2 /s 3 ) versus time at 40 m 
AGL from the 0000 UTC initial data. 
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TKE error 


TKE Error vs. Time 
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Figure 12 - Simulated averaged turbulence kinetic energy error (m 2 /s 2 ) versus time at 40 
m AGL from the 1200 UTC initial data. 
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TKE Error vs. Time 
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Figure 13 - Simulated averaged turbulence kinetic energy error (m 2 /s 2 ) versus time at 40 
m AGL from the 0000 UTC initial data. 
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